Revealing the Supernova— Gamma-Ray Burst Connection with TeV Neutrinos 
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Gamma-ray bursts (GRBs) are rare but powerful explosions displaying highly relativistic jets. It 
has been suggested that a significant fraction of the much more frequent core-collapse supernovae 
are accompanied by comparably energetic but mildly relativistic jets, which would indicate an 
underlying supernova-GRB connection. We calculate the neutrino spectra from the decays of pions 
and kaons produced in jets in supernovae, and show that the kaon contribution is dominant and 
provides a sharp break near 20 TeV, which is a sensitive probe of the conditions inside the jet. For a 
supernova at 10 Mpc, 30 events above 100 GeV are expected in a 10 s burst in the IceCube detector. 



PACS numbers: 96.40.Tv, 97.60.Bw, 98.70.Sa 

Long duration gamma-ray bursts (GRBs) have been 
found to be tightly connected with core-collapse super- 
novae [jj] (short duration GRBs may be caused by an- 
other mechanism such as compact star mergers 0). Al- 
though many mysteries remain, this strongly indicates 
that the central remnants of the core-collapse event, most 
likely black holes, drive the observed relativistic (bulk 
Lorentz factor Tf, > 100) jets of GRBs. High-energy neu- 
trinos from accelerated protons in GRB internal shocks 
are predicted to be detectable at future 1 km 3 Ccrcnkov 

r-i 

detectors such as IceCube |3( ■ (If the protons escape and 
collide with external material, both the delayed gamma 
rays from neutral pions and the neutrinos from charged 
pions could be observable The jet signature, but 

more mildly relativistic, may be common to supernovae, 
which are much more frequent than GRBs, even correct- 
ing for the effects of jet opening angle on observability. If 
a significant fraction of core-collapse supernovae are ac- 
companied by jets with ~ 3, perhaps ~ 1% according 
to late-time radio observations [5j, then neutrinos could 
be the only prompt signature of these hidden sources. 

In addition to being more frequent, mildly relativis- 
tic jets are expected to be much more baryon-rich (a 
"dirty fireball" ) . Both properties work quite positively 
for neutrino detectability. Recently, Razzaque, Meszaros, 
and Waxman (RMW) developed a model of high-energy 
(TeV) neutrino emission from jets with Tf, ~ 3 Q. Be- 
cause of the low Lorentz factor and high baryon den- 
sity, collisions among accelerated protons (pp) occur ef- 
ficiently, making pions that decay into neutrinos, and a 
nearby core-collapse supernova at 3 Mpc is predicted to 
be detectable at IceCube 0. If those high-energy neu- 
trinos are detected, and correlated with an optical super- 
nova, then it would strongly and directly show the pres- 
ence of a mildly relativistic jet with significant kinetic en- 
ergy, as well as potentially being the first detection of ex- 
tragalactic neutrinos. Data from many supernovae would 
give the distribution of the jet Lorentz factor, providing 
important insight into the supernova-GRB connection. 



We extend the RMW model and significantly improve 
the detection prospects; most importantly, we consider 
the kaon contribution, since compared to pions, kaons 
have several advantages. Due to their larger mass and 
shorter lifetime, kaons undergo less energy loss before 
decaying into neutrinos. This, and the relatively larger 
energy transferred to the daughter neutrinos, means that 
the neutrino spectrum from kaons is harder and more de- 
tectable than that from pions. The neutrino spectrum 
from kaon decay has a sharp drop near 20 TeV, which is 
a sensitive diagnostic of the acceleration mechanism and 
the conditions inside the jet. Since many detected events 
are expected (e.g., ~ 30 at 10 Mpc in IceCube), our re- 
sults greatly extend the range of detectability, which in- 
creases the frequency as ~ d 3 , where d is the distance 
of the furthest detectable objects. While these models 
have many uncertainties, our results significantly improve 
their detectability, such that even the existing AMANDA 
detector can provide important constraints. 

Jet Dynamics. — We briefly summarize the RMW 
model, using their notation The jet kinetic energy 
is set to be Ej — 3 x 10 51 erg, which is typical for GRBs. 
Since we discuss a baryon-rich jet, we use the mildly rel- 
ativistic value for the bulk Lorentz factor Tf, = 3, and 
an assumed opening angle of 9j ~ T^ 1 = 0.3. By anal- 
ogy with observed GRBs, it is natural to set the vari- 
ability timescale of the central object as t v = 0.1 s. 
The internal shocks due to shell collisions then occur 
at a radius rj — 2T%ct v — 5 x 10 10 cm, smaller than 
a typical stellar radius. The comoving number densities 
of electrons and protons inside the jet are n' e = n' p = 
Ejj( y 2'KQ]r]Ylm v c z tj) = 4 x 10 20 cm" 3 , where tj = 10 s 
is the typical GRB jet duration; the superscript ' repre- 
sents quantities in the comoving frame of the jet. It is 
assumed that fractions e e = es = 0.1 of the jet kinetic 
energy are converted into relativistic electrons and mag- 
netic fields, also by analogy to GRBs. These electrons 
lose energy immediately by synchrotron radiation. The 
synchrotron photons, however, thermalize because of the 
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high optical depth, n' e aT^bct v = 2 x 10 6 , where <jt is the 
Thomson cross section, which makes the jets invisible 
with prompt gamma rays, unlike GRBs. As a result, the 
magnetic field strength, photon temperature and number 
density are given by B' = [AesEj / '(9]r]Tlctj)} x / 2 = 10 9 
G, El, = [15(hc) :i eeE : j/{2ir 3 9]r]r 2 b ct : j)} 1 / 4 = 4 keV, and 
nl, = E j e e /(2ne]r]Tlct j El 1 ) = 8 x 10 24 cm" 3 . It is as- 
sumed that the internal shocks accelerate protons with 
a spectrum ~ E~ 2 , normalized to the jet total energy. 
The maximum proton energy is obtained by compar- 
ing the proton acceleration timescale t' &cc ~ E' p /eB' = 
1CP 12 s (E' p /1 GeV) |7| with the energy-loss timescales. 
RMW assume that radiative cooling by the synchrotron 
and Bethe-Heitler (pry — * pe + e~) processes dominates, 
and that the maximum proton energy is 2 x 10 6 GeV. 

Proton and Meson Cooling. — We describe our exten- 
sion of the RMW model, including the maximum accel- 
eration and meson cooling arguments. At energies be- 
low the photopion production (hereafter pj) threshold 
E' pi th = 0.3/El, GeV 2 = 7 x 10 4 GeV, the proton ac- 
celeration timescale is much shorter than any energy-loss 
timescale. Above the P7 threshold, we find that the most 
competitive cooling mechanism is the jyy process itself, 
due to the very high photon density. Assuming 15% en- 
ergy lost from an incident proton in each pry interaction, 
AE' p = 0.l5E' p , and using a pl = 5 x KT 28 cm 2 ||, we 
obtain a cooling timescale t' pl = E' p j (c.Op 1 n' 1 AE' p ) = 6 x 
1(T 8 s. Equating t' acc = t' pj , we obtain E p = 5x 10 4 GeV, 
slightly less than the threshold energy of the pj interac- 
tion. Thus when the p^j interaction becomes accessible, it 
prevents further acceleration, so we take E' p max = 7 x 10 4 
GeV, in contrast to RMW. As discussed below, we focus 
on the break in the kaon-decay neutrino spectrum as a 
direct observable of the maximum proton energy. 

Accelerated protons produce mesons efficiently via 
pp interactions, since the optical depth is r' = 
n' p a pp T b ct v — 2 x 10 5 , where a pp — 5 x 1CT 26 cm 2 ||. 
The meson multiplicity in each pp interaction is taken 
to be 1 for pions, and 0.1 for kaons; this is the right 
ratio 0, but we have been conservative in the normal- 
ization to focus on the most energetic mesons. We as- 
sume that the mesons are produced with 20% of the 
parent proton energy, so that they follow the original 
spectrum of accelerated protons with power-law index 
—2. They cool, however, by radiative (synchrotron ra- 
diation and inverse Compton scattering off thermal pho- 
tons) and hadronic (irp and Kp) processes. The cool- 
ing timescales are t' rc = 3m 4 c 3 /[4<7tto 2 -E' '({/' + U' B )] and 
t' hc = E' /{c(Jhn' p AE') for radiative and hadronic cool- 
ing, where m and E' are the meson (jr or K) mass and 
energy, — E'n' and U' B = B' 2 /8it are the energy 
densities of photon and magnetic fields, ah is the cross 
section for meson-proton collisions, and AE' is the en- 
ergy lost by the incident meson in each collision. Adopt- 
ing AE' = 0.8E' 10] and a h = 5 x 10~ 26 cm 2 for 
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FIG. 1: Diagram for the neutrino spectrum from pion or kaon 
decays. The spectral break energies are given by Eqs. Q 
and for pions, and by Eqs. and Q for kaons. These 
spectral features are smeared out in detection (see Fig. 



pions and kaons t' hc is energy independent, while 
t' rc oc E'^ 1 . Since the total cooling timescale is given by 
t'^ 1 = t' h ~ c + 1'^ 1 , it is dominated by hadronic cooling at 
lower energies, and radiative cooling at higher energies. 

If mesons decay faster than they cool, then the daugh- 
ter neutrinos maintain the spectrum shape; otherwise, 
the spectrum becomes steeper. Below the cooling break 
energy E'^\ obtained by equating j't = t' c (^ t' hc ), 
where 7' and r are the meson Lorentz factor and proper 
lifetime, there is no suppression because mesons imme- 
diately decay after production. Above E'^\ the cool- 
ing suppression factor is the ratio of two timescales, i.e., 
t'cll' T - When the energy-loss process is dominated by 
hadronic cooling, the suppression factor is oc E'~ l . When 
it is dominated by radiative cooling, above an energy 
E^ b , it is oc E'~ 2 ; we obtain E^ b by solving t' hc = t' rc . 
RMW neglected hadronic cooling of mesons; it is not a 
large effect for pions, but it is important for kaons. 

Neutrino Spectrum. — Charged pions and kaons, which 
are assumed to carry 20% of the original proton energy, 



decay into neutrinos through ■n :il ,K ± 



with branching ratios 100% and 63%. The neutrino en- 
ergy in the observer frame is related to the parent meson 
energy in the jet rest frame as follows: E v = T^E'^/A 
and E v = TbE' K /2. The Lorentz factor represents the 
Doppler boosting effect, and 4 and 2 in the denomina- 
tors reflect the fraction of the parent energy conveyed 
by the daughter neutrino in the case of pion and kaon 
decays. Secondary neutrinos from muon decays are irrel- 
evant since muons immediately undergo radiative cool- 
ing |6j. The energies and densities here are similar to 
the case of neutrino production in Earth's atmosphere, 
which makes the calculation more robust. 

Figure n is a diagram for the neutrino spectrum from 
meson decays, characterized by two spectral breaks and 
a maximum neutrino energy. For neutrinos from pion 
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decay, we obtain the cooling break energies: 

ETj$ = 30 GeV, (1) 
= 100 GeV, (2) 

corresponding to E'^ b and E'^ b . The dependence on the 
jet parameters is given by Ej T b 9Hjt 2 and (e e + eB) _1 r b 
for Eqs. ||TJ and (J2J, respectively. We note that the first 
break energy is strongly sensitive to the value of r& (it 
is less severe if one assumes Tf, ~ 9j , following RMW). 
This means that the model is quite uncertain, but at the 
same time, that the detection of neutrinos could precisely 
constrain the Lorentz factor of the jet. 

The neutrino spectrum from kaon decays is much more 
favorable, for three reasons. First, radiative cooling is 
much less efficient than for pions, since kaons are heavier 
and the radiative cooling timescale is t' rc oc m 4 . Sec- 
ond, the kaon lifetime is a factor ~ 2 shorter. Third, a 
larger mass also shortens the particle lifetime because of 
a smaller Lorentz factor at fixed energy. Thus the cooling 
breaks of kaons occur at much higher energies: 

= 200 GeV, (3) 

E *6b = 20 ' 000 GeV > ( 4 ) 

where the scaling is the same as Eqs. and (jSJ. The 
maximum energy £„, max = ^bE' K ma ^/2 is only slightly 
above the second break for a canonical parameter set, al- 
though this could be changed for other parameter choices. 
Measurement of the sharp edge of the neutrino spectrum 
would be a sensitive test of the maximum proton energy, 
and hence the physical conditions in the jet. 

Neutrino Burst Detection. — We first estimate the nor- 
malization of the neutrino spectrum, evaluating the flu- 
ence at the first break energy, F u> q = F v (E^ cb ). Assum- 
ing efficient energy conversion from protons to mesons, 
and that half of the mesons arc charged, we obtain 
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where d is the source distance, (n) is the meson multiplic- 
ity (1 for pions and 0.1 for kaons), B v is the branching 
ratio of the decay into neutrino mode (1 for pions and 0.6 
for kaons), and the factor ln(E' p>max / E' in ) normalizes 
the proton spectrum to the jet energy. For canonical pa- 
rameter choices and for a nearby source at d = 10 Mpc, 
F„, becomes 5 x 10~ 2 and 5 x 10~ 5 GeV -1 cm~ 2 , for 
neutrinos from pion and kaon decays, respectively. The 
parameter dependence is E^T b 14 9~ 6 t~ 2 t~ i d^ 2 . 

We calculated the expected signal from one supernova 
neutrino burst, using the code ANIS (All Neutrino In- 
teraction Generator) 0]. We neglect the effects of neu- 
trino oscillations, as they are below the uncertainties of 
the model. Figure Pfa) shows the event spectrum from 
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FIG. 2: (a) Event spectrum of neutrino-produced muons from 
a supernova at 10 Mpc in a 1 km 3 detector. Contributions 
from and K decays are shown as dotted and dashed 
curves, and the total as a solid curve. The atmospheric neu- 
trino background is shown for comparison; it is evaluated for 
1 day and within a circle of 3° radius, (b) The same, but 
cumulative event number above a given energy. 

the muon neutrinos and antineutrinos from a supernova 
at 10 Mpc, and in Fig. |2fb), we show the yields above 
a given energy. We used a detector effective area of 1 
km 2 , which is reasonable for IceCube in the case of up- 
going muons 0. We took into account the muon range, 
which effectively enlarges the detector volume, and eval- 
uated the muon energy when it enters the detector if it is 
produced outside, or at the production point otherwise. 
Since the spectrum of neutrinos from pions falls steeply, 
their expected event spectrum is also steep, and there- 
fore, if we lower the threshold, many more events would 
be expected, as shown in Fig. I3b). The spectrum of 
neutrinos from kaons, on the other hand, is much flatter, 
making them the dominant component at high energies. 

If we take 100 GeV as the threshold muon energy for 
IceCube (for a transient point source), we expect about 
30 events from a core-collapse supernova at 10 Mpc, 
mostly from kaons. (If the proton spectral index is not 
—2.0, but is instead —1.5 or —2.5, the expected number 
of events is 40 or 3, respectively.) These events cluster 
in a 10 s time bin and a ~ 3° angular bin surrounding 
the supernova, which allows very strong rejection of at- 
mospheric neutrino backgrounds. If the source is farther 
and the expected number only a few events, then we may 
use a more conservative time bin, e.g., a 1-day bin cor- 
related with optical observations, considering the time 
uncertainty between the neutrino burst and an optical 
supernova. The atmospheric neutrino background for 1 
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day and within a circle of 3° radius is also shown in Fig. 
|5J The spectral break corresponding to the maximum 
proton energy might be detectable around ~ 20 TeV, 
if the source is very close or the jet energy Ej larger than 
assumed. Recalling that the proton acceleration is lim- 
ited by the fry interaction, this probes the photon num- 
ber density in the jet, an essential quantity of the model. 
If the neutrinos come from below the detector, it may 
be possible to reduce the detector threshold, so that the 
large yield due to pion decays could be detectable too. 

Discussion and Conclusions. — The RMW model, 
while speculative, is an specific and intriguing proposal 
for a common thread connecting GRBs and core-collapse 
supernovae, namely the presence of jets with energies 
around 3 x 10 51 erg. For GRBs, which are rare events, the 
jets would be highly relativistic and revealed; for super- 
novae, which are frequent, the jets would be mildly rela- 
tivistic and hidden, except from neutrino telescopes. The 
most basic features of the RMW model are a proton spec- 
trum falling as E~ 2 , carrying most of the jet energy, and 
a suitable target density for neutrino production. In anal- 
ogy to the observed properties of GRBs, and evidence for 
a supernova-GRB connection, these requirements do not 
seem unreasonable. What is significant is that these con- 
siderations can be directly and easily constrained with 
neutrino detectors, lessening the dependence on theory. 

The key remaining issue is the fraction of supernovae 
which have these jets. Some type Ic supernovae do sug- 
gest the presence of mildly relativistic jets, based on their 
late-time radio emission; the fraction of all core-collapse 
supernovae with jets is perhaps ~ 1% 0, significantly 
larger than the fraction of supernovae with highly rela- 
tivistic jets, i.e., GRBs. It is possible that mildly rel- 
ativistic jets are in fact much more common, but are 
completely choked in the hydrogen envelopes of type II 
supernovae (RMW assumed that nearly all supernovae 
have jets). If such mildly relativistic jets accompanied SN 
1987A and were directed toward Earth, the Kamiokande- 
II and 1MB detectors would have seen many events, the 
small detector sizes being more than compensated by the 
closer distance. This suggests that such jets did not exist 
in SN 1987A; additionally, it was dark in the radio [l^ . 

In any case, the nearby core-collapse supernova rate is 
high enough to allow testing these models soon, especially 
if type II supernovae have jets which can only be revealed 
by neutrinos. Within 10 Mpc, the rate of core-collapse 
supernovae is more than 1 yr _1 , with a large contribu- 
tion from galaxies around 3-4 Mpc [x3 | . At 3 Mpc, we 
expect about 300 events in a 10 s burst in IceCube, which 
would be very dramatic. The electron neutrino shower 
channel, where the detected energy more faithfully re- 
veals the neutrino energy, may be especially helpful |l5j . 
Even in the already-operational AMANDA, though the 
effective area is smaller than that of IceCube, we expect 
~ 10 events. Since these events would arrive within 10 
s, from a specific point source, the background rejection 



should be excellent, and AMANDA could already con- 
strain these models. 

At larger distances, the galaxy distribution is smooth 
enough to allow simple scaling: the number of sources 
increases as d 3 , while the neutrino signal per supernova 
decreases as d~ 2 . For example, at 20 Mpc, the expected 
number of neutrino events in IceCube is still several, and 
the total supernova rate is greater than about 10 yr _1 . 
That distance contains the Virgo cluster (which increases 
the supernova rate somewhat beyond the simple scaling) 
in the Northern hemisphere, which is an attractive target 
for South Pole neutrino telescopes. The effect of beaming 
is to reduce the source frequency at a given distance, but 
to increase the source fluence in a compensating way. 
Given the large assumed opening angle, the probability 
of having the jets directed towards Earth is relatively 
large, around ~ 10%. Even taking into account that 
perhaps not all supernovae have jets, this model predicts 
rich detection prospects for IceCube. The detection of a 
prompt burst of high-energy neutrinos would reveal the 
time and direction of the core collapse event, which would 
be useful for forecasting the optical supernova, and for 
searching for a gravitational wave signal in coincidence. 
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